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Abstract 

A domain in a Langmuir monolayer can be expected to have a shape that 
reflects the textural anisotropy of the material it contains. This paper ex- 
plores the consequences of XY-like ordering. It is found that an extension of 
the Wulff construction allows for the calculation of two-dimensional domain 
shapes when each segment of the perimeter has an energy that depends both 
on its orientation and its location. This generalized Wulff construction and 
newly-derived exact expressions for the order parameter texture in a circular 
domain lead to results for the shape of large domains. The most striking 
result is that, under general conditions, such domains will inevitably develop 
a cusp. We show that the development of cusps is mathematically related to 
phase transitions. The present approach is equivalent to a Landau mean-field 
version of the theory. 

PACS numbers: 61.30. Cz, 68.10.-m, 68.35.Md, 82.65.Dp 
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I. INTRODUCTION 



The regularities in the shapes of crystalline materials have long been understood to be 
macroscopic expressions of the positional order of crystals at the atomic level. The angles 
between the faces of a crystal are, for instance, structural invariants dependent only on a 
certain set of integers characterizing the faces (Miller indices [I]). In a classic paper 0, 
Wulff developed a geometrical construction allowing for the determination of the shape of a 
crystal, provided one knows the values of the surface energy of the crystal as one varies the 
Miller indices. One "draws" planes with every possible set of Miller indices. The distance of 
each plane from a fixed point in space is proportional to the energy per unit area of a surface 
parallel to that plane. The inner envelope of the planes is then the equilibrium shape of a 
finite piece of crystal. 

Key to the Wulff construction is the dependence of the interfacial energy on orientation. 
Anisotropic surface energies are, however, not restricted to the crystalline solids. A portion 
of the surface of a liquid-crystalline mesophase has, in general, a surface energy whose value 
depends on its orientation relative to the optical axis. The principal difference between 
this system and a crystal is that liquid crystals may have crystallographic directions that 
are easily deformable (e.g. smectics) or they may lack crystallographic order altogether — 
and only exhibit orientational order (e.g. nematics). It seems obvious that the internal 
softness of liquid crystals and like systems will suppress such characteristic surface features 
of solid crystals as facets and edges. In the case of an extremely soft system the optical 
axis ought to adjust itself so as to minimize the effective surface energy for all orientations. 
A sample is then expected to relax to a spherical shape to minimize its overall interfacial 
energy. According to this argument, one should expect to encounter edge-like features only 
in liquid-crystalline materials that are relatively rigid. As we will see in this paper, this last 
conclusion is, in fact, quite incorrect. Edges can be seen even in very soft liquid crystals. 

The possibility of sharp edges, or cusps, in liquid crystalline materials was first suggested 
by Herring in the early 1950's. He utilized the Wulff construction to demonstrate that liquid 
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crystals with rigid orientational order can have interfaces with sharp cusps, provided that the 
surface energy is sufficiently anisotropic. The resulting cusp angle are no longer geometrical 
invariants but instead depend on material parameters. It was also noted by Burton, Cabrera 
and Frank (BCF) in 1951 that the Wulff construction acquires important simplifications in 
two dimensions (D = 2) J/J. Its derivation becomes a simple exercise in the calculus of 
variations, in contrast to the situation in three dimensions, for which there is no entirely 
rigorous and general demonstration of the correctness of the Wulff construction ||. To 
implement the Wulff construction in two dimensions, one solves a simple linear differential 
equation, while in three dimensions the only known method is the geometrical construction 
described above. 

Anisotropic shapes of samples in thermal equilibrium are, indeed, encountered in exper- 
iments on D = 2 mesophases. The systems in question are Langmuir monolayers consisting 



of surfactant molecules at an air- water interface ||11|| . The phase diagram of such materials 
generally contain, in addition to the D = 2 analogue of the solid liquid and gaseous phases, 
the so-called liquid condensed (LC) phase, which is liquid-crystalline. The LC phases are 
two dimensional anisotropic liquids in which the direction of the tilted hydrophobic tails 
define an anisotropy axis, c. Because of the anisotropy between the polar head groups and 
the non-polar tail there is no c — > —c symmetry (as is the case for nematic liquid crystals). 

At the phase boundary between the LC and isotropic liquid phases, stable coexistence 
droplets are observed with shapes that are, in general, non-circular. Recently, polarized 
fluorescence microscopy studies of coexistence droplets of pentadecanoic acid with a linear 
dimension in excess of 25 /im have reported that the boundaries of these droplets contain a 



cusp JT2|. Such a cusp in illustrated in Fig. la. The interior angle of the cusp was found to 
increase with the size of the domain. 

The existence of cusps on the boundary of these droplets would not seem to be too 
surprising, as droplets in surfactant systems have been found to contain toplogical defects 
in the pattern — or texture — of the anisotropy axis. A defect located at the boundary could 
well produce a cusp. However, Brewster angle microscopy studies — which allow for detailed 



visualization of textures — on these materials have revealed defect-free texures of the "virtual 
boojum" type (see Fig. lb) with the cusp in the sample shape located opposite the 
focus of the boojum. 

The determination of the shape of these systems demands a simultaneous minimization 
of the free energy with respect to the shape and the internal structure of the sample. In 
other words, we must find the analogue of the Wulff construction for deformable media. 
The question of finding such a generalized Wulff construction for Langmuir layers belongs 
to a wider class of shape minimization problems, as encountered, for instance, in ferrofluid 
droplets 0, where the shape energy depends on the internal structure of the sample far from 
the surface. Analytical treatment of such problems beyond perturbation theory has proven 
difficult because the coupling between different parts of the surface, as mediated by the 
deformable bulk structure, is highly non-local. As suggested by BCF [0] , shape calculations 
are more easily carried out in D = 2. We will demonstrate in this article that a generalized 
Wulff construction can be found in D = 2 for the case when the internal structure of the 
sample is describable by an X — Y model. As discussed below, Langmuir monolayers in the 
LC have an order parameter which, under certain conditions, reduces to an X — Y model. 
We will show that cusp singularities ought to be generic features of 2D drops with an X — Y 
order parameter. 

We will start in Section |J with a discussion of the free energy minimization for a fixed 
sample shape (a circle) and the resulting textures. In Section || we will keep the texture fixed 
and minimize the free energy with respect to sampl shape, following the method of BCF 
0. Simultaneous minimization of the free energy with respect to both shape and sample 
texture can be done perturbatively, as discussed in Section |II[ Our general method, based 
on complex function theory, is discussed in Section [VT]. A brief summary of the method has 
been published elsewhere [||. 
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II. FREE ENERGY OF XY DROPLETS 



In this Section we define the Hamiltonian and examine the order parameter textures. We 
will assume a rigid circular boundary, so free energy minimization simply involves finding 
the lowest energy texture consistent with the boundary energy. The order parameter is 
taken to be the unit vector c = (cos 0, sin 0) with an associated effective Hamiltonian 

R[e(x,y)} = f ^Ve 2 dxdy+f a{6-Q)dS. (2.1) 

J interior 2 J boundary 

The first term in Eq. ( |2.1| ) describes the free energy cost associated wih a posistion depen- 
dence of the order parameter c(r). The coefficient k is, for the XY model, the so-called 
spin-wave stiffness. For Langmuir monolayers in the LC phase, k corresponds to the Frank 
constant for the special case that the so-called "splay" and "bend" Frank constants are 
equal. The second term on the right hand side or Eq. (|2.1| ) represents the "surface" energy 
of the domain. The variable 9 is the angle that the unit normal makes with respect to the 
axis to which the director angle is referred. 

If the boundary has a fixed shape, minimization of the energy in Eq. (|2.1|) implies the 
following two extremum equations 

V 2 0(x,i/) = O, (2.2) 

and 

on 



Eq. ( [2.2|) applies in the bulk, and Eq. fl2.3|) at the surface of the domain. The derivative 
dQ/dn in Eq. (|2.3| ) is along the surface normal, and cr'{x) = da(x)/dx. The bulk extremum 
equation, Eq. (|2.2| ), requires a Q(x,y) that is a harmonic function. The general solution of 
Eq. ( |2.2|) in D = 2 can then be written, for the case at hand, as: 

0(x, y) = - (f(x + ly) - f(x - iy)) , (2.4) 
% 
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with f(z) an arbitrary analytic function of z = x + iy. As an example of this method, 
expand f(z) in a Taylor series: 



n=l 



(2.5) 



Then, using Eq. Q2.4p, we find 



0(x, y) = aij/ + a 2 ya; + a 3 ?/(3x 2 - y 2 ) + 
1 

z 



(2.6a) 
(2.6b) 



which is the most general analytic texture with c(r = 0) = x and ^(x, y) = — c y (x, —y). 
The most general form for the surface energy is 

oo 

(7(9 -S) = 5>„cos(n(#-0)). (2.7) 

n=0 

Physically, the coefficient a can be identified as the isotropic line tension. The coefficient 
ai measures the lowest order surface anisotropy. If a\ < 0, then the associated surface 
energy favors a c-vector along the outward boundary normal, while for a\ > it favors a c 
lying along the inward normal. For Langmuir monolayers, a\ is in general nonzero, but for 
D = 2 nematic liquid crystals a\ = by symmetry. In that case we must go to the next 
term (n = 2). If a 2 < then the surface energy favors a normal orientation for c at the 
boundary without distinguishong whether it is outward or inward, while a 2 > favors a 
parallel orientation for c. 

Under the assumption of a circular domain, the boundary condition, Eq. (|2.3|) , becomes 



k [e ie f(e w ) - e- i9 f(e- ie ) 

" e »9 e n(-/(e' s )+/(e-")) _ e -« e n(/(e»)-/( e - i( )) 



D oo 

+ -7T 1^ na n 

Z n=0 



0. 



(2.8) 



Eq. ( |2.8| ) can be solved by iteration when a n <C k/R . For instance, if only a and a\ are 
non-zero, then we find a Taylor series for f(z): 



fi \ aiR ° 1 ( aiR ° V , r\( 3\ 



(2.9) 



The choice of the function f(z) is thus imposed by the solution of Eq. fl2.3|) . Comparing 
Eqs. ( p.9| ) and ( |2.6[ ) one can verify ct\ = —aiR /2n and a 2 = — (ciiR q /2k) 2 . The resulting 
texture is shown in Fig. lb. It has a mathematical singularity outside the sample, which is 
called a "virtual" boojum. Textures of this type are, indeed, frequently encountered in both 
Langmuir monolayers and on the surface of smectic C* liquid crystals. The distance of this 
virtual singularity in the textural structure from the center of the domain, Rb, is obtained 
by extrapolating lines of constant to their point of intersection. Inserting Eq. ( p.9|) into 
Eq.( [2.4|) and performing the required extrapolation, one obtains 

R B = 2K/ ai R . (2.10) 

Note that as ai — > 0, the boojum recedes to infinity and the texture becomes uniform 

For the case n = 2, the 2D nematic drop, we must set a\ = 0. If only a 2 is non-zero, 

then interative solution of Eq. ( |2.8|) yields 

„. . 1 a 2 z 2 
f(z) ^ — , 

so «i = and a 2 = —a 2 /RoK. The resulting texture now has two virtual singularities, but 
this time they are located along the y axis. 

The above result for a\ and a 2 , for the case a\ ^ 0, in the perturbation region a\ -C k/Ro 
suggests a Taylor expansion of f(z) with the coefficients a n = - (aiR /2n) n . This is indeed 
an exact solution of Eq. ( |2.8| ). In fact, whenever the boundary energy has the special form 
o~(9 — 0) = a n cosn(6 — 0), i.e. whenever only one a n is nonzero (besides ao), an exact 
solution of Eq. ( |2.8[ ) can be found for f(z): 

f(z) = -\og(l-a n z n ), (2.11) 
n 



nan Ro/k , 
a n Rl = , " 0/ _ = . (2.12) 



with 



1 + y 1 + (na n R / k) 

A demonstration that Eqs. ( |2.11| ) and ( |2.12| ) indeed yield a solution to the boundary con 
dition Eq. (|2.8|) is contained in Appendix A. 



If n = 1, then the exact solution corresponds to a virtual boojum, lying a distance Rb 
away from the center of the domain, where 



As the ratio aiRo/n goes to zero, the boojum again retreats to infinity As cliRo/k goes to 
infinity, corresponding to a very strong anisotropic surface energy, or a very large domain, 
Rb — > Ro + n/a%. The boojum thus approaches the sample in this regime, but the spacing 
remains finite in the limit Rq — > oo. 

When n is greater than one, the exact solution for the texture is equivalent to the 
texture produced by n singularities lying outside of the domain (see Figs. 2b-d). The 
virtual singularities are no longer boojums. The singularities for n > 2 are "fractionally 
charged" in the sense that we do not recover the starting orientation of the director field if 
we perform a circuit around the singularity. For n > 2 the resulting singularity is equivalent 
to a nematic disclination. There are singular lines attached to the singularities where 
jumps by 2ir(n — 2)/n, but these lines do not intersect the sample. 

Textures corresponding to a%, . . . , 04 ^ are displayed in Fig. 2. For clarity of presen- 
tation, the singularities are placed on the boundary of the domains. 



We now turn to the second part of the problem: minimizing the free energy with respect 
to sample shape for a rigid texture. We will first review the analysis of BCF of the 
D = 2 Wulff construction, and then we will use their method to see under what condtitions 
we ought to expect to encounter sample shapes with cusps. The texture is assumed to be 
uniform in the analysis immediately following. 

We start by introducing a parameterization of a two dimensional curve. This parameter- 
ization expresses the co-ordinates of a point on the curve in terms of the angle between the 
tangent to the curve and the x-axis. If this angle is 9 and the distance between the tangent 
line and the origin is R(9), then 





III. WULFF CONSTRUCTION IN TWO DIMENSIONS 



S 



x{9) = R{9) cos(#) - sin(0) 
y (0) = i?(0)sin(0)+^cos(0). 



(3.1a) 
(3.1b) 



The quantities 9 and R(9) are displayed in Fig. 3. It is relatively straightforward to express 
9 and R{9) in terms of x and y. One has 



cot(0) = 

ax 



R{9) 

1 + {dy/dx 

The following relations follow immediately from Eqs. ( |3~T 



x 



dx 



(3.2) 



dx 
d9 



sin(0) 



d# 2 



d y (q\ 

- = cos(0) 



R{9) + 



d 2 R(9) 
d9 2 



(3.3a) 
(3.3b) 



and 



dS_ 
~d9 



( dx 



+ 



d9 



R{9) + 



d 2 R(9) 



d9 2 



(3.4) 



Furthermore, the area inside a closed curve is simply expressed as an integral over 9: 

d 2 R{9) 



A = ~ j [xdy - ydx] = -j R{9) 



R(9) + 



d9 2 



(3.5) 



The integral in Eq. ( p.5|) is taken counterclockwise around the curve. 

Suppose, now, that one wishes to minimize the boundary energy of a two dimensional 
domain having an anisotropic surface energy, u{9). Here, the angle 9 is both the angle that 
the boundary's unit normal makes with respect to the x-axis and the angle that parame- 
terizes the bounding curve in the parameterization of Eq. (|3.1|) . This minimization is to be 
achieved subject to the constraint that the total enclosed area is a constant. Using Lagrange 
multipliers, we arrive at the following extremum equation: 
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5R(9) 



R(9') 



d 2 R{9') 



d 2 a{9) 
d9 2 



R{9) + 



d9> 2 
d 2 R[9) 
d9 2 



- \m 



R(9') + 



d 2 R{9') 
d9> 2 



d9' 



(3.6) 



or, 



R{9) + 



d 2 R{8) 
d9 2 



1 

A 



a(6) + 



d 2 a(9) 
d9 2 



(3.7) 



The solution of the above equation is 



1 



R{0) = -a(9) + d cos(6) + C 2 sin(0). 
A 



(3.8) 



According to Eq. ( |3.8| ), apart from the C\ and C2 terms, the minimum energy shape has 
a bounding curve such that R(9) is proportional to the anisotropic surface tension, a (9). 
Since R{9) is the distance from the origin to the tangent of the bounding curve, we have 
recovered the Wulff construction precisely. It can be shown, by direct substitution into Eqs. 
fl3.1|) , that the only effect of the additive sine and cosine terms is to translate the domain 
without changing its shape. 



As a first example of the use of Eqs. (3.1) and (|3.8|) , consider again the case in which 
only ao and a\ are finite in the fourier expansion, Eq. fl2.7|) , of <y{9). Setting = in Eq. 
then gives 

(3.9) 



a(9) = a + a,\ cos(0). 

Inserting Eq. ( |3.9| ) into the right hand side of Eq. 

d 2 R{9) 1 " 



we obtain 



R(9) + 



d9 2 



A 
1 
A 



ao + a% cos 9 + -r^(a + a± cos 1 
d9 2 



a . 



(3.10) 



The solution of Eq. (|3.10|) is, then, identical to that for the case a± = 0, which is a circle. 
According to Eq. (|3.10|) a simple cos 9 contribution to the anisotropic surface energy leads 
to no distortion of the shape of the domain. In fact, we show in Appendix B that in general 
a contribution to a(9) that is proportional to cos 9 cannot affect the sample shape for rigid 
textures and, in particular, cannot produce any cusps. 
Next, consider a surface energy of the form: 



a (9) = a Q e^ cos{e) . 
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(3.11) 



Using Eq. (|3.11|) and inserting the result in Eq. ( |3.8|) with C\ = C2 = gives 



x (e) = \cos{o) + p sm(ey\ e ^ cosW , (3.12a) 

y (0) = [1 - /3cos(0)] sm{e)e^ cos ^. (3.12b) 

Fig. 4 displays the boundaries for (3 = 0.5, 1 and 1.5. For (3 = 1.5, the bounding curve has 
a "swallowtail" feature. According to the rule that one keeps only the inner envelope, the 
proper prescription is to amputate this feature and retain only the left-hand portion of the 
curve, which now has a cusp. The development of swallowtail singularities in a family of 
one-parameter curves is a familiar feature of the theory of catastrophes. It is, for instance, 
encountered in a set of curves evolving according to the Huyghens construction of physical 
optics, where swallowtail singularities are associated with the development of caustics. In 
our case, the curve parameter is /3, and this parameter must have a finite value for swallowtail 
singularities to develop. In Appendix C we prove that the swallowtail feature appears first 
at 13 = 1, so the sample boudary is smooth for (3 < 1. 

The key difference between the two anisotropic boundary energies, Eq. ( |3.10| ) and Eq. 
( p.ll| ), is that the latter boundary energy contains higher order harmonics. Indeed, if we 
include the higher harmonics of Eq. 



a (9) = J2 a nCos(n9). (3.13) 



n=0 



and insert Eq. ( [3.13j ) into Eq. ( |3.7|) 



R{9) + = £ a n (l - n 2 ) cos(n0). (3.14) 

then, as soon as Fourier components with n > 2 are included in Eq. ( |3.14| ) we find a 
distorted boundary, but as long as the a n coefficients with n > 2 remain small compared 
to do, we find no cusp singularity. The swallowtail features cannot in general be obtained 
perturbatively. In other words, the anisotropic contribution to a (8) must be comparable to 
the isotropic line energy ao before any singularities in the bounding shape can develop. This 
would suggest that the experimentally observed cusp feature can only be explained if one 
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assumes a highly anisotropic line tension. In the Section |V| we will use perturbation theory 
to show why this conclusion is quite incorrect. 



IV. CONNECTION WITH PHASE TRANSITIONS 



The example discussed in the previous Section also serves to illustrate the mathemat- 
ical connection between the onset of a cusp and a thermodynamic phase transition. This 
connection is more readily developed if one recasts the problem of calculating the minimum 
energy bounding surface into more conventional notation. In this notation the procedure 
giving rise to a cusp produces expressions identical to those encountered in the standard </> 4 
model of a symmetry-breaking phase transition. 

We begin by writing the expression for the domain's bounding curve in the form x = f(y)- 
This single- valued function describes the domain in the immediate vicinity of y = 0, where 
the cusp will occur in the case at hand. The orient at ion- dependent surface tension has the 
form 



The surface tension of the portion of the boundary curve lying in the top half of the x — y 
plane is 



Taking the functional derivative of the above expression with respect to f(y), we obtain the 
following extremum equation, wich applies at the end-point y = 0. 



a(9) 



e (3/y/l+(dx/dy) 2 



(4.1) 




(4.2) 




(4.3b) 



(4.3a) 
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When (3 < 1 the only solution to the above equation is /'(0) = 0, while when (3 > 1 there 



are the additional solutions corresponding to yl + f'(0) 2 = (3, or, returning to angular 
variables, (3 = 1/ cos(# cusp ). 

To gain further insight into the mathematical nature of the onset of the cusp, we expand 
the term in brackets on the right hand side of Eq. fl4.3aQ . If f'(0) is small and (3 ~ 1, then 



F(f(0)) = y/l + fWet'^WW 



1 + (1 - (3)f'(0) 2 + V (0) 4 + O (f(0) 6 , ((3 - l)f(0) 



(4.4) 



The quantity F (/'(0)) is the surface energy per unit length along the y-axis, at the point the 



cusp develops. The right hand side of Eq. (4.4) has precisely the form of a Ginzburg-Landau 



like mean field theory. The combination 1 — (3 plays the role of the reduced temperature 
and /'(0) is the order parameter. The expression possesses a local minimum at /'(0) = 
when (3 < 1, which becomes a local maximum as (3 passes through one. The onset of the 
cusp corresonds to the system's choosing on the the non-zero values of /'(0) associated with 
the local minima at 



/' (0) = ±2 J (3 - 1 (4.5) 



Only one of the solutions of Eq. (|4.5|) is physically relevant. According to the Wulff con- 
struction, the sign of /'(0) must be the one associated with an outward-pointing cusp. This 
implies a small symmetry-breaking term in the theory, at this point unidentified in our 
analysis. 

In general, if the surface tension cr(8), has a maximum at 6 = 0, then Eq. ( f4.4|) can be 
written as 

F (/(0)) = & [(7(0) + <r"(0)] (/'(0)) 2 + C 2 (/'(0)) 4 . (4.6) 

Cusps appear when cr(0)+cr"(0) changes sign. This is consistent with the results of Appendix 
C, where it is demonstrated that cusps appear when R(0) + R"(0) changes sign. 
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V. SOFT LIMIT: EFFECTIVE SURFACE TENSION 



Before developing the general formalism, we will first consider the limiting case of very 
small k, for which the problem simplifies significantly. In this limit we can define an effective 
surface tension which can be used in the two-dimensional Wulff construction, as applied to 
rigid materials. In other words, we include the deformability of the material by a re-definition 
of the surface energy. 

In the soft limit the texture of the sample deforms itself to respond to the surface energy 
anisotropy. We start by setting k = 0. The texture can then adjust itself freely to minimize 
the anisotropic line tension. Let <To be this minimum value. Since the line tension is a 
constant, and since the textural energy is zero, the sample shape must be circular. The 
line tension is minimized when the director c lies along the outward-directed normal to the 
circle, so c is in the radial direction along the circle perimeter. The associated texture must 
be a solution of Eq. ( |2.2|) obeying this boundary condition. In general, the solution will 
have one or more singularities. We will focus on the solution of Laplace's equation with a 
singularity on the circle perimeter: the boojum texture: 



Now, let K be small but finite. Two things must happen: (i) the c-director exerts a 
torque on the boundary which is then deformed away from a perfectly circular shape, and 
(ii) the c-director is no longer perfectly along the boundary normal. We will first keep the 
texture fixed at the at Eq. (|5.1| ), allowing the sample shape to relax, and then we will allow 
the texture to relax and reconsider the sample shape. 

First, use Gauss's law to rewrite the total energy as a line integral over the boundary 
assuming 6 to be a solution of Eq. (|2.2|) : 




(5.1) 




(5.2) 



The quantity cr e ff given by 
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(5.3) 



then appears as an effective line tension. The function 0(s) is, however, a functional of 
the boundary shape so the right hand side of Eq. ( |5.3| ) is really a non-local expression and 
cannot be simply interpreted as a line tension. If, nevertheless, we use use Eq. Q5.1|) in Eq. 
r3f) we find: 



K sin 9 

&efi = 7-^°-, n + °"0- 



(5.4) 



2R 1-cosfl 

The resulting effective surface tension is an analytic, but aperiodic function of 9. The line 
tension anisotropy only depends on the dimensionless parameter V = n/aoR, with R the 
sample radius. It can be inserted into the Wulff construction and leads to a cusp at 9 = 0, 
with a cusp angle proportional to T. 

We now redo this calculation while allowing for textural relaxation. First, assume a 
circular boundary. Allow G to deviate from the boojum texture, and expand the line energy 
around its minimum a^: 



6(x, y) = Q B (x, y) + w(x, y) 
a{9-Q) ^a + ^(9-<d) 2 . 



(5.5) 
(5.6) 



To find w we must solve V 2 u> = subject to the boundary conditions Eq. ( |2.3| ). In terms 
of w: 



K 



dw 
Or 



+ a"(0)w(s) = -k 



dS, 



r=R 



dr 



(5.7) 



r=R 



The solution of this mixed boundary condition problem for w can be written as a Fourier 
expansion in 9: 



w(r, 



2J a m r m smm9 



(5.8) 



m=l 



with coefficients: 



7 ra"(0)R m {l + ^ ) ) J-« 
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dr 



(5.9) 



r=R 



Inserting Eqs. (|5.8|) and ( p^9|) into Eq. (|5T2|) yields: 



F = R 



d9\a - -a^6w(r = R,9) + 0(w A ] 



(5.10) 



where we have used Eq. ( |5.7| ). We neglect the fourth and higher order terms in w. Inserting 
the solution Eq. ( |5.8| ) into the line integral expression for F gives, as expected, a non-local 
expression: 

,a de B (6') 



F = R I d9{a + 9 I d9'K(9, 9' 
with a kernel K(9,9'): 



dr 



(5.11) 



r=RJ 



Ts/n n'\ K sin m.6 1 sin m.6 1 ' 
m9) = 2^^ 



i _i_ rriK 



(5.12) 



In the large R limit, this kernel reduces to two delta functions, at 9 = ±9', while for finite 
R, these two delta functions broaden by an amount 59 of order k/<j"R. To second order in 
k/R we can neglect the spreading of the delta functions. In the large R limit, the non-local 
line energy thus becomes a local line tension. Using Eq. Q5.12| ) in ( |5.11| ), one obtains: 



F = R 
R 



2 or 



(5.13) 



OG 



Surprisingly, we recover the "naive" line tension of Eq. (|5.4j ). If we now allow the shape 
to relax, we find the shape discussed earlier. This shape relaxation must then be used in a 
recalculation of w, but in the large R limit such corrections are higher order in 1/R. We 
conclude that in the limit of small k and large R the sample shape is expected to have a 
cusp with an excluded angle proportional to k/ctq-R. 



VI. GENERALIZATION OF THE WULFF CONSTRUCTION 



We now relax the condition of small k and develop our method, which allows for a 
determination of the sample shape, even when the a e s of Eq. (p.3|) is truly non-local. Our 



16 



method will be a generalization of the BCF procedure and the results of Section 0. We 
will assume that the texture always obeys Laplace's equation, Eq. ( |2.2| ), with the director 
angle specified by a complex function f(z) through Eq. (|2.4|). The function f(z) is, for a 
given sample shape, determined by Eq. ( |2.8| ). The remaining problem is now to minimize the 
sample free energy with respect to sample shape. We will use the BCF parameterization R(9) 
for the sample shape. To find the optimal sample shape, we must compute the variational 
derivative of both surface and textural energy and equate it to the variational derivative of 
the sample area A with respect to R{9) times the Lagrange multiplier A — in direct analogy 
with Eq. ( |3.6|) . We will, in the following, always assume near-circular shapes. 

We start with the variational derivative of the textural energy. Imagine a two dimensional 
domain containing an order parameter described by the director angle Q(x,y). If Q(x,y) is 
given by Eq. (|2.4|) then 

V6| 2 = 4f'{x + iy)f'(x-iy), (6.1) 
and the bulk contribution to the energy in Eq. ( [2.1Q is, then, 

2k J F(x,y(x))dx, (6.2) 
where y(x) denotes the boundary line energy, and where 

F(x, y) = [ V f\x + iy')f{x - iy') dy'. (6.3) 



The textural energy, Eq. ( |6.2|) is explicitly dependent on the shape of the boundary, and, 
again, can be treated as a non-local contribution to the line tension. Using Eqs. ( j3.3|) and 
( |3.4| ) we rewrite the infinitesimal dx as — sin^rfs, with s(9) the arclength. The derivative 
with respect to R{9) of the textural energy Eq. (|6~^ ) is then 

d 2k [ F{x {6'), y {6')) sin 6' ds{6') 



5R{6) 

= 2k f (x (9) + iy (9)) f (x (9) - iy (9)) . (6.4) 

We now turn to the variational derivative of the surface energy, Eq. ( ^.7| ). Using the 
Fourier expansion, Eq. (|2.7| ), for a{9 — 6) we find: 
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j a(e-e) ds 



5R(6) 

Re < V] a n e~ ind e n ^ x+iy ^~ n ^ x ~ iv ^ 



n(n + l)(x + iy)f'(x + zy) + n(n — l)(x — iy)f{x — iy) + 1 — n 2 j . 

(6.5) 



The algebraic steps in the derivation of Eqs. ( |6.4|) and ( |6.5D are outlined in Appendix D. 

The two variational derivatives appear to be forbiddingly intricate. However, our complex 
function f(z) is not just any function; it must obey Eq. Q2.8|) , and we use this to simplify 
Eqs. ( |6.4| ) and ( |6.5| ). First, we rescale lengths in the circular domain so that its radius is 
unity, as described in Section |I|. On the boundary of the domain, one can then replace 
x + iy, respectively exp(m#) by z, respectively z n (where \z\ = 1) and the quantities x — iy, 
respectively exp(— inO) by 1/z, respectively l/z n . Equation ( |2.8| ) which determines f(z) 
then simplifies to: 



K 



z"e 



n(-f(z)+f(l/z)) _ 



z -n e n(f(z)-f(l/z)) 



0. 



(6.6) 



Dividing both sides by z and integrating yields 



K 



f 



x z \x 



dx 



E 



+ K 



z n e -nf{z)+nf(l/z) _|_ z ~n e nf(z)-nf(l/ z) 

xf(x)~-f (-)} \f(x) + - 2 f (-) 
x \x/i L x z \x J 



dx (6.7) 



Using Eqs. ( |6.6| ) and ( |6.7| ) in Eqs. ( |6.4| ) and ( |6.5| ) we find that the functional derivative 
with respect to R(0) of the energy of the domain, as given by the sum of the right hand 
sides of Eqs. fl6.4p and ( |6.5p , is equal to J-(z), with 



T{z) 



ZK 



dz 



^w-^(i)] + k W'wi 2+k [^'(i) 



-K 



m + fl- 



wf{wf-—f(- 



dw, 



(6.8) 
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up to an unimportant constant. 

The equation determining the optimal shape is now found by repeating the steps between 
Eqs. ( |3.5| ) and (|3.6|), with the result: 

m + = X [a ° Ro + m] ' (6 - 9) 

where <Jq is the isotropic contribution to the surface tension, and where F{9) is given by Eq. 
Q6.8D with z replaced by e td . 

Equations (|6.8| ) and (|6.9| ), which hold as long as the domain boundary is not too far 
deformed from a circular shape, form the basis of our method. The sample shape is deter- 
mined in three steps: (i) solve Eq. ( |6.6|) to obtain f(z) for a given anisotropic line tension 
a(9 — Q), (ii) compute J 7 (9) from Eq. ( |6.8|) to obtain the effective line tension, and (iii) solve 
the BCF formula, Eq. (|6.9| ) with JF \9) to find R{9). We will now consider some examples 
of the application of our method. 



A. = cfq + a\ cos 4> 

Our first example is the case in which a n = for n > 2 in the Fourier expansion, Eq. 
( p.7|) . Recall that for an imposed, rigid, uniform texture the minimum energy shape was 
a perfect circle (Eq. ( |3. 10|) ) , and that for an imposed circular domain shape the minimum 



energy texture was the virtual boojum (Eq. ( [2 . 1 1|) with n — 1). To implement our recipe 
for finding the shape which minimizes the total free energy, we first note that the virtual 
boojum texture: 

f(z) = log(l - a x z) (6.10) 



is an exact solution of Eq. ( |6.6j ), with a\ given by Eq. ( p,12| ). If we then insert f(z) into 



J-(z) (Eq. (|6.8j)) we find T(z) = constant. The resulting shape equation 

R{9) + = T [*oRo + const.] (6.11) 
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is the equation for a perfect circle. We conclude that the virtual boojum texture in a circular 
domain is a free energy minimum, both with respect to domain shape and domain texture. 
We have yet to find a simple physical argument that explains why the perfect circle remains 
a free energy minimum for this anisotropic form of the line tension. 



B. (j{4>) = cq + a n cos n<j) {n > 1) 



Let us now try the same procedure for n > 1. First, assume a perfect circular sample. 

The associated texture then follows from f(z) = - log(l — a n z n ), which is a solution of Eq. 

If we now compute T{z) and insert the result in the shape equation we find 
d 2 R(6) 



R(6) + 



d6< 



0Q 



(n 



1 — n 



n 



K 

+ — 



1 - a n e in6 
1 

+ 



[1 - a n e 
1 — n 



ind\2 



n 



1 - a n e- ine 



a„e 



-in0\2 



(6.12) 



The circular shape is no longer a free energy minimum. In the limit of very small domains, 
with a R /n <C 1, Eq. (|6.11|) reduces to 



R(9) + 



d 2 R(6) R 



d6 2 



(TO 



a + (1 — n 2 )a n cos nO 



(6.13) 



(6.14) 



with solution 

R(0) ^ R (l + —cosne) . 

\ CO J 

The domain shape has an n'th harmonic shape deformation on top of the circular shape. 
Note that the domain shape is independent of the Frank constant k. This result is, of course, 
not quite exact, since f(z) will, for a non-circular shape, no longer be given by Eq. fl2.ll ), 
but the resulting corrections are higher order in a n /o~o and o~qRq/ 'k. 



C. a(4>) = cfo + a\ cos (j) + (Z2 cos 2(j> 

We have seen under A that for a pure cos0 line tension anisotropy we can find a free 
energy minimum in the perfectly circular shape. Since we have an exact solution, we can use 
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perturbation theory to asses the effect of higher order harmonics. We thus include one more 
term — n = 2, with a 2 <C a\ — and recompute the shape. This particular form of the surface 
term has been argued to be relevant for Langmuir films |I6| . From our earlier results one 
would naively expect a nearly circular sample shape with a small n = 2 correction. However, 
things turn out a little differently 

If we take a 2 <C a\ then it will be possible to expand f(z) about its a 2 = form. Writing 



f{z) =\og(l - a z) + f 1 (z), 



(6.15) 



with a = aiRo given by Eq. ( |2.12|) , the boundary condition Eq. ( |6.6| ) becomes, to first 
order in a 2 and f±, 

1 — az~ 
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Ro 



zf[{z) - -f[ (- 
z \z 



z — a 
1 — az 



+ 



z — a 



^ (I 

= 0. 



z — a 



+ 



1 — az' 


a-2 
2 




' z — a ' 
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'1 — az' 
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z — a _ 




.1 — az_ 




. z — a _ 





(6.16) 



We separate Eq. ( |6.16| ) into two equations: 



Oi 

2 



— a 



a 



z — a 



+ 



1 — az' 


do 


' z — a ' 


z — a _ 


+ i 


.1 — az. 



0. 



(6.17) 



.1 — az 

and an identical equation with z replaced by 1/z. Each of those equations is a simple linear 
first order differential equation, solvable by standard methods. One finds 



a 2 a z — a 



V a /q-° a ) zt ~ a dt. 



(6.18) 



a\ 1 — a 2 1 — az J 1 — azt 

Computing T{z) with the use of f(z) = log(l — az) + fi(z) one finds the shape equation 



1 K 

r[z) = 



f l (z) + 2a(l-az) J — 



-dx 



az) 



(6.19) 



where R(z) = r(z) + l/r(z). The final reduction of the result for R(9) consists of algebraic 
manipulations. First, we note that the quantity a = a\Ro as given by Eq. ( |2.12j ) approaches 
1 as Rq — > oo. This means that, in the limit that Rq is large we can replace 1 — a by 
k/qiRq. Furthermore, inspection reveals that in the limit k/cliRo 1 and a 2 a\ the first 
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order contribution to x(6),y(6) is negligible compared to the zeroth order one except in the 
immediate vicinity of 9 = 0. After a series of changes of variables that take this behavior 
into account we arrive at the following expression for R(6) 

Ik di. 



R{6) 



-Re 



(Xo-Ro o,\ 

where the dimensionless function X(^) is given by 



(6.20) 



1 + ii) 



1 — iip Jo 1 — iip + y 



dy 



-4(1 -iip) 



oo e -y 
oo p~ y 



log(l — iip + y) — log(l — iip) 



+ 



I -iip 2(1 



+2(1 - #) / — - log(l - # + y) + log(l - #) + 

jo y 1 — i 



dy. 



(6.21) 



y" { i — iip 

A graph of the real part of the function X(^>) is displayed in Fig. 5. The following 
properties of the function play an important role in the behavior of the boundary 

1. Re [X(?/>)] is an even, nonsingular function of its argument, and can thus be expanded 
in an power series in ip. 



2. Re [X(-0)] has a negative second derivative at ip — 0. 

3. Re [X(^)] possesses a minumum at ip ~ 2.74. 



1. Onset of a cusp 



As discussed in Section |TV] and Appendix C, the onset of a cusp is signaled by a change 
in sign of the function R(8) + d 2 R(8)/d8 2 . From the previous arguments it follows that this 
is most likely to happen at 9 = 0. Making use of our result, Eq. Q6.20 ), for R(6) we find 



R(0) + 



d 2 R{6) 



de 2 



1 + 



Ik a 2 



ReX(O) + 



1d20,\R[ 



X"(0). 



(6.22) 



=0 c"o-Ro a i °~0 K 

In the limit of vary small a>2 the expression on the right hand side will be positive unless 

the domain radius Rq exceeds a threshold value, the leading contribution to which is given 

by 
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— 2a ia2 |X" 

<7nK 

= 0.356—, (6.23) 

where the number in the last equality above follows from a numerical evaluation of the 
second derivative. When the sample radius is slightly in excess of -Rthrcshoid the excluded 
cusp angle, ip (the difference between the cusp angle and 2n), obeys 



tjj OC JR - ^threshold- (6.24) 



As in the case of the example discussed in Section [TV], the behavior of the cusp angle at 
onset obeys a power law consistent with a <f) 4 mean field theory. A graph of the excluded 
angle of the cusp as a function of domain radius can be found in Fig. 6. The domain radius 
is expressed in units of k/cli. In that case, the only other adjustable parameter is the ratio 
a-ijoQ. In the Figure, that ratio is equal to 0.05. 

2. Large R limit 

The behavior of the excluded angle when R is very large can also be determined by 
inspection of Eq. ( |6.20 ). As the symmetry of the domain is the same in this case as in the 



illustrative case discussed in Section [TTT] and in Appendix C, we locate a cusp by searching 
for a solution to the equation y(6) = 0, with 9 small but not equal to zero. Substituting the 
right hand side of Eq. ( |6.2U| ) into Eq. ( |3.1b| ) we find that, for small 9 and large Rq, with 



9R finite, 

= Za^Ro g j Re[^'(9R a 1 /K)] } ^ 
cr K \ 8RqCLi/k, J 

Because of the properties of the function X(-?/>) enumerated above, the ratio in curly 
brackets in Eq. ( |6.25| ) is well-behaved as a function of the variable OR^ai/ ' k. Neglected in 



Eq. (|6.25| ) are terms of zeroth order in Rq . 

According to Eq. flOg ), y(6) = when either 9 = or Re [X' (ORo^/k)] = 0. Finite 



cusp angles correspond to the second case. Because Re [X(^)] when ip = 2.74, we conclude 
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that 9 = 2.74 (n/R Q ai). The interior angle of the cusp is then equal to 7r — 5A8k/ciiR (see 
Appendix C) and the excluded angle of the cusp, defined as 7r minus the interior angle, is 
equal to 5A8k/ciiRo. The excluded angle decays as 1/Rq and is independent of 0,2 when Rq 
is sufficiently large. 

Fig. 7 shows what the domain boundary looks like when a cusp is induced by the 
02 cos 20 term. The energy parameters have been set so that the inequalities relied upon 
in this Section are satisfied. In particular, g^/oo = 0.25 and Rqcii/k = 5. Note that a 
swallowtail appears, just as in the illustrative example cited previously. 



VII. CONCLUSIONS AND COMPARISON WITH EXPERIMENT 

In summary, we have demonstrated that a two-dimensional domain containing a de- 
formable medium describable by an X — Y model has a purely circular boundary, as long 
as the surface anisotropy term is of the form a\ cos 0, where is the angle between the 
X — Y vector and the unit normal to the bounding curve. This remarkable result no longer 
holds if one introduces the small anisotropic surface energy 02 cos 20. In certain regimes the 
surface deformation is singular; a cusp appears when the sample size exceeds a threshold 
determined by the coefficients of the anisotropic and the bulk Frank constant, k. The cusp 
angle ultimately decays with increasing sample size as one over the domain radius. The 
threshold domain radius is proportional to <7oft/ a i a 2j where <7o is the isotropic bounary en- 
ergy. Note that the more deformable the material in the domain is (i.e. the smaller the 
value of k) the more likely it is to develop a cusp — which contradicts naive intuition based 
on the appearance of facets in crystalline materials. 

How does our theory compare with the experimentally measured cusp angles? In Fig. 8 



we show a set of cusp angles measured as a function of sample size by Schwartz, et al [12 
for pentadecanoic acid. The data has been fitted to the curve displayed in Figure 6. The 
ratio a^jo has been set equal to 0.05. No attempt was made to adjust that ratio so as to 
optimize the fit, but the quatities a± and k were effectively varied by adjusting the vertical 
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and horizontal scales. Note that the measured cusp angles show no evidence of a sharp onset 
at a threshold radius. 

An even more serious discrepancy between the data and our results lies in the fact that 
the best fit to the asymptotic power law decay of the measured cusp angle with domain size 
is cusp oc -Ro°' 3 . Our mean field theory gives a decay as 1/Rq. One possible origin of the 
discrepancy is that we did not identify the correct operator as the one that "breaks" the 
degeneracy of the pure ai cos ip boundary energy. Langmuir monolayers have Frank constants 
K\ and K3 for splay and bend that are not identical [0|. We have assumed in Eq. ( p . 1\) 
that Ki = K3. The bulk equation satisfied by the director angle O is now nonlinear. Is the 
solution to this equation related in any way to the virtual boojum solution that plays such 
an important role in the determination of the domain's shape? Preliminary work indicates 
that cusps also appear if we let K\ 7^ k 3 . Another possibility is that thermal fluctuations 



lead to a renormalization of the effective anisotropic boundary energies |1y| . This possibility 
is under current study. 

The cusp singularities discussed in this paper should be generic features of boundary lines 
in Langmuir layers and two-dimensional liquid crystalline materials in general. Indeed, defect 
lines in hexatic Langmuir monolayers frequently exhibit a scalloped appearance (JTjJ, which 
may have the same origin as the cusps discussed here. The extension of the work reported 
here to more general sample shapes is not straightforward, with regard to the family of 
exact solutions for the texture we have found, an obvious question is whether some variation 
on a conformal transformation allows one to generate from it the texture appropriate to a 
non-circular domain. It is readily verified that a simple conformal transformation will not 
simultaneously deform the domain's bounding curve and appropriately alter the texture. 
This is because of the nonlinear nature of the boundary conditions. However, given that we 
have obtained a set of exact solutions for the texture in the presence of realistic boundary 
conditions, an extension of the applicability of those solutions would be highly desirable. 
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APPENDIX A: THE VIRTUAL BOOJUM SOLUTION 



In this Appendix we show that Eqs. ( 2.11|) and ( p. 12 ) represent exact solutions to the 
boundary value problem represented by Eq. ( |2.8| ), with one nonzero a n . The demonstration 
is by direct substitution. We have 



and 
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-ind 
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ind 
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-in9 



1 — a n e 



-ind 



+ 



1 QLiyy^t- 1 Oif-f^G 



(Al) 



(A2) 



According to Eq. (|2.8| ), the sum of the left hand sides of Eqs. (|A1|) and (|A2|) equals zero. 
The two right hand sides sum to zero if 



1 2k 
a n — = 0. 

a n a n Ro 



(A3) 



The solution to this equation for the parameter a n that has the proper limiting behavior is 
displayed in Eq. (|2.12|) . 
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APPENDIX B: STABILITY OF THE BOUNDARY AGAINST a{9) oc cos(0) 



Consider an anisotropic surface tension having the form a (9) = Acos(9). The total 
boundary energy of a closed boundary of arbitrary shape is 

j> c ■ hds. (Bl) 

The quantity c in the above expression is a constant vector, h is the unit normal to the 
boundary and ds is the infinitesimal element of length along the boundary. In two dimensions 
we can write 

c • hds — z ■ (c x ds) , (B2) 

where z is the unit vector out of the plane and ds is the directed infinitesimal length element 
along the boundary This means that the boundary energy is 

(B3) 

The closed integral on the right hand side of Eq. ( |B3| ) is always equal to zero. This means 
that a a(6) oc cos(#) has absolutely no effect on the shape of a two dimensional domain. 



ds) x c 



APPENDIX C: ON THE APPEARANCE OF CUSPS 

To demonstrate exactly how the implementation of the Wulff construction in two dimen- 
sions leads to cusps on boundaries we consider the equation for the bounding curve generated 
by Eqs. ( |3.12| ) — that is, the boundary when the anisotropic surface tension is given by Eq. 



( [3. 11| ). Because of the symmetry of the domain, as illustrated in Figs 4, a cusp lying on the 
rightmost edge of the domain represents a solution to the equation y{9) = 0, subject to the 
additional condition 9^0. According to Eq. (|3.12b|) , the above requirements are met when 



P cos9 = 1. (CI) 

This equality can only be satisfied if (3 > 1. By simple geometry, we find that the interior 
angle of the cusp is equal to ir — 29 c , where 9 C = arccos(l//3) is the angle satisfying Eq. (|Cl| ). 
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Fig. 4c depicts the boundary when (3 = 1.5, while Fig. 4b show how the domain looks 
when f3 is equal to its critical value of 1. 

Of interest is the "swallowtail" appendage attached to the cusp in the former case. 
The construction of the inner envelope corresponding to the equilibrium boundary shape is 
completed by the amputation of this appendage. Although the behavior of x{9) and y(6) in 
the tail is, thus, nominally irrelevant to the shape of the domain, a brief discussion of the 
properties of x and y in the tail region yields useful information regarding the mathematical 
signals that accompany the onset of a cusp. Consider, for instance, the cusps in the tail, at 
either end of the "trailing" edge. Using Eqs. (|3.3j ), we find dy/dx = — cot(6 ) ). Thus, there 
will be no discontinuity in the slope of the curve given by Eqs. ( |3.1| ) as long as 9 varies 
continuously. On the other hand, if the combination R{9) + d 2 R(6)/d6 2 changes sign, then 
both x(6) and y(8) reverse directions with increasing 6, leading to a cusp at which both 
segments of the curve meet tangentially. Such cusps represent the only type of discontinuity 
that can appear in the boundary generated by the construction outlined in Eqs. ( |3.1| ) — in 
the absence of an amputation at a point at which the boundary crosses itself. If we assume 
that any cusp generated by the Wulff construction will make its appearance accompanied by 
a swallowtail, then the onset of a cusp is signaled by a change in sign of R{9) + d 2 R(9)/d9 2 . 
This sign change first appears at the point on the boundary at which the cusp develops. 



In the implementation of the generalized version of the Wulff construction one is, char- 
acteristically, faced with the task of taking the variational derivative with respect to R{9) 
of an integral of the form 



where A (sin 6, cos 9) is a polynomial in sin 8 and cos 8, f (x (8) ,y (6)) is a polynomial in 
x{6) and y{6) and dS = [R{0) + d 2 R(9)/d9 2 } d6 is the infinitesimal length element along the 



APPENDIX D: ON VARIATIONAL DERIVATIVES 




(Dl) 
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boundary. While the relatively simple dependence of x and y on 9 allows the variational 
derivatives to be implemented in a straightforward manner, the complete and accurate 
determination of the derivative with respect to R{9) of an integral in which the functions 
f(x, y) and A(sin 9, cos 9) have any but the simplest form becomes extremely tedious, in the 
absence of a stratagem that allows for the evaluation of functional derivatives of non-trivial 
integrands. Fortunately, such a stratagem exists. 

Consider the following version of the integral above: 



sin 9f(x, y) dS. 



(D2) 



First, note that 



sin 9 dS = sin#— 
d9 



sin# 



R(9) + 



d 2 R{9) 
d9 2 



dx{9) 
d9 ' 



(D3) 



The integral can thus be rewritten as 



- I f(x(0),v(0))%dO. 



(D4) 



Now, we take the derivative 8/8R{9) of the above integral: 
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(D5) 



Now, given Eqs. (|3.1| ), we have 

5x 5 



8(9 - 9') cos 9 -5' {9 -9') sin^, 



(D6) 
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and 



5y 



8R{6') 



6(9-9') sm9 + 8'(9-9')cos9. 



(D7) 



With the use of the above two relations and Eqs. ( |3.1| ) the right hand side of Eq. QD5I 
becomes 
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Thus, 



8R(9) 



f df 
J sm9f(x,y)dS = — 



R(6) + 



d 2 R{9) 
d9 2 



(D9) 



A similar set of manipulations yields the result 

5 r . .„ df 



8R{9) 



J cos9f(x,y) dS = — 



R{9) + 



d 2 R(9) 
d9 2 



(D10) 



Note that the final results hold no matter what form is taken by the function f(x, y). The 
method described above can be fruitfully applied when the integrands are more complicated 
than the one above. 

It is important to note that the steps leading to the final result, Eq. ( P10|) , included 
integrations by parts, in which the "perfect derivative" terms were neglected. If the domain 
has a boundary that is free of singularities — especially cusps — one can easily justify this. 
However, when cusps are present, so that periodic boundary conditions on an integral around 
the circumference of a domain wall cannot be assumed, more attention to those terms is 
called for. The absence of a detailed analysis of the effects of a cusp in the boundary on the 
results of an integration by parts is a gap in the development in this work. Such an analysis 
is clearly called for. 
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FIGURES 

FIG. 1. A drawing of the kind of cusped domains seen in experiments. Such domains contain 
liquid condensed regions in an environment of the liquid expanded phase of a Langmuir monolayer. 
Figure la depicts the domain as it is seen in fluorescence microscopy. Figure lb adds the "virtual 
boojum" texture that is conjectured to be responsible for the cusp. 

FIG. 2. The various virtual singularity textures that emerge from exact solutions to Eqs. (2.2) 
and (2.3) when only one of the a n 's in the surface energy in non-zero and the domain is circular. 
Figure 2a: n = 1, 2b: n = 2, 2c: n = 3, 2d: n = 4. The singularities sit outside the domain, except 
in the limit of inifinite k. For ease of depiction, they are shown on the domain boundary. 

FIG. 3. The quantities 9 and R(9), as defined in Eqs. (3.1). 

FIG. 4. The shape of a domain, when the surface tension has the anisotropic dependence on 
9, e P co < e \ and x{9) and y{9) are given by Eqs. (3.12). As indicated in the Figure, the three cases 
depicted are f3 = 0.5, = 1 and (3 = 1.5. 

FIG. 5. The function Re [X(^)], as defined in Eq.(6.21) 

FIG. 6. A plot of the excluded angle = tt — ip- m versus the radius of the domain. The 
radius of the domain is expressed in units of K/a±. The ratio 02/170 has been set equal to 0.05. 
The quantities are defined in Section VI. The quantity ipi n is the interior angle of the cusp, and is 
shown in Figure 7 

FIG. 7. The domain when 0,2/0-0 = 0.25 and RqQi/k = 5. For definitions of the quantities, see 
Section VI. The interior angle, ipin, is indicated in the Figure. 
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FIG. 8. A comparison between the curve for the excluded cusp angle displayed in Figure 6 
and data obtained by Schwartz, Tsao and Knobler (see ref. 12). The axes are both logarithmic 
and scales are adjusted to obtain a by-eye fit. No attempt was made to optimize agreement by 
adjusting energy parameters in the theoretical result 
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